Fast growth at low temperature in vacancy-mediated phase-separation 
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We study the phase-separation dynamics of a two- 
dimensional Ising model where A and B particles can only 
exchange position with a vacancy. In a wide range of temper- 
atures the kinetics is dominated, during a long preasymptotic 
regime, by diffusion processes of particles along domain inter- 
faces. The dynamical exponent z associated to this mecha- 
nism differs from the one usually expected for Kawasaki dy- 
namics and is shown to assume different values depending 
on temperature and relative AB concentration. At low tem- 
peratures, in particular, domains grow as t 1/2 , for equal AB 
volume fractions. 



When binary alloys are quenched into the unstable 
region of the phase-diagram, the two species A and B 
segregate into coarsening domains of typical size L(t). 
The simplest model for the description of the dynami- 
cal behavior of these systems is the Ising model evolving 
with Kawasaki A-B exchange rule [1|. Despite a great 
deal of work, mostly numerical [||-3|, the understand- 
ing of this phase-separation kinetics is still not com- 
pletely satisfactory. For example the expected asymp- 
totic growth L(t) ~ t x l z with z = 3 due to the 
evaporation-condensation mechanism (bulk diffusion) |t]] 
has not been directly observed in simulations so far, be- 
cause of extremely long transients. This poor compre- 
hension is also in contrast with the wide recent interest 
in submonolayer ordering in epitaxial systems, for which 
the two-dimensional Ising model is the simplest descrip- 
tion [||. In real alloys direct A-B exchange is extremely 
unlikely; swaps between particles and vacancies are in- 
stead the fundamental mechanism leading towards equi- 
libration p0| . Investigations of this kind of dynamics 
in the Ising model have been performed pl|^2| showing, 
at intermediate temperatures, the convergence towards a 
late stage z = 3 exponent to be enhanced with respect to 
the Kawasaki case. However, despite these equal asymp- 
totic exponents, the issue of the universality of large scale 
properties for the two kinds of dynamics is far from be- 
ing clarified due to a poor understanding of the basic 
processes governing growth. Several mechanisms have 
been proposed f|,[l^] for both cases, predicting a wealth 
of preasymptotic exponents, but no clear identification 



of them has been found in previous numerical studies. 

In this Rapid Communication we investigate vacancy- 
mediated dynamics in a wide range of temperatures and 
relative concentrations. We identify long preasymptotic 
regimes during which the growth exponent has well de- 
fined values different from the value 1/3 expected for the 
evaporation-condensation mechanism. In particular, at 
low temperatures, the growth exponent is larger than at 
intermediate temperatures and, for equal volume frac- 
tion of the two species, the remarkably fast growth law 
L(t) ~ i 1 / 2 is found. This rich phenomenology, observed 
numerically, is interpreted in a unified and consistent 
framework by means of simple scaling arguments. A 
prominent role is played by diffusion processes along the 
edges of segregating domains. This mechanism couples 
to the rough or faceted nature of the interfaces at dif- 
ferent temperatures and to the domain morphologies at 
different concentrations, to give rise to a variety of growth 
laws. The fact that particle-exchange is mediated by va- 
cancy motion turns out also to be crucial, producing a 
notably different preasymptotic pattern with respect to 
the Kawasaki dynamics. 

We consider an Ising-type model with Hamiltonian 
H = —JJ2<ij> s i s j where Si — ±1,0 corresponds to 
A or B particles and holes, respectively, and the inter- 
action is between nearest neighbors. The only micro- 
scopic motion allowed is the exchange between a vacancy 
and one of the nearest neighbor particles. At time zero 
the system is instantaneously brought from an initially 
high-temperature uncorrelated state to the final temper- 
ature T < T c . We perform Monte Carlo simulations of 
the model evolution with the Metropolis algorithm on a 
square lattice of size 512 x 512 with one single vacancy. 
This hole density is in the range of typical values for 
real alloys |l3|. Time is measured as the number of at- 
tempted vacancy steps. We calculate the characteristic 
domain size L(t) as the position of the first zero of the 
correlation function averaged along the vertical and hor- 
izontal axes. Other measures of the domain size, as the 
inverse of the first moment of the radially averaged struc- 
ture factor, give similar results. The behavior of L(t) for 
symmetric quenches is shown in Fig. [|J here and in the 
following temperatures are rescaled by the critical tem- 
perature T c — 2.269 . . . J/fcs fl4| ■ For very short times 
the vacancy has not yet visited all parts of the lattice. 
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This causes the unusual behavior of L(t) for t < 10 8 . For 
longer times the evolution varies remarkably depending 
on temperature. At T = 0.6, the domain size quickly 
starts growing as t 1 ^ 3 |l2|]. At low temperature (T = 0.2) 
the effective exponent becomes larger than 1/3 and set- 
tles to 1/2 during the last decade considered. For in- 
termediate temperatures an even more surprising result 
occurs. The value of L(t) for T = 0.45 becomes, af- 
ter about 10 9 steps, larger than the same quantity for 
T = 0.6. On much longer times, the same occurs also for 
T = 0.3. This result is at odds with the common wisdom 
according to which the higher the temperature, the faster 
the separation process. 

For asymmetric quenches, such that the fraction of A 
particles is 80%, results are illustrated in Fig. Again 
for low temperature the effective growth exponent is 
larger than for higher T in the range of times consid- 
ered. In particular we find that at T — 0.6, 1/z *» 0.27 
while, at T = 0.3, 1/z is very close to 1/3. 

In the following we present scaling arguments explain- 
ing the observed behaviors in a unified framework. The 
basic ingredient is that, in the time range considered, 
diffusion of particles along domain interfaces is the rel- 
evant process. This physical mechanism alone allows a 
consistent interpretation of all the phenomenology. Our 
strategy to determine the dynamical exponent z is the 
following: We consider first a droplet of A particles of 
typical size L deformed from its equilibrium shape. By 
denoting with t eq the time needed to relax the perturba- 
tion we show that t eq ~ L@ and determine (5 for the dif- 
ferent temperatures considered in the simulations. Then 
we deduce a relationship between f3 and z both for sym- 
metric and asymmetric compositions. 

Relaxation time of a droplet deformation. The 
exponent depends on the roughness of the interfaces 
which, in turn, is determined by temperature and the 
scale considered. The roughening temperature for the 
d = 2 Ising model is zero and therefore interfaces are 
always rough on large scales. However, domain walls ap- 
pear faceted at finite T on length scales smaller than a 
typical length L C (T) which can be estimated to be of 
order exp(2J/T)/2 jl5|] . For sufficiently low tempera- 
tures, then, L(t) < L C (T) in an appreciable time inter- 
val and interfaces are straight on the characteristic scale 
L(t). This is exactly what occurs in our simulations at 
T = 0.2. For such a temperature, L C (T = 0.2) w 41 and, 
for all times considered in Fig. |l| L(t) < L C (T). When 
T = 0.6, instead, L C (T) s=s 2 and in the phase-separation 
process domain walls are always rough. We consider the 
case with flat interfaces first. The basic process is the 
following (Fig. ||): A vacancy is captured on the outer 
part of the droplet interface and wanders along it until 
it switches to the inner part by moving through a kink 
located at P. A passage through another kink in Q and a 
successive evaporation bring the hole in the bulk of the B 
phase again. The net effect is the transfer of an A parti- 



cle from Q to P. Since the vacancy moves without energy 
costs along the facets, the duration of the process is deter- 
mined by the activated events (evaporation and passage 
through kinks) which take a time r ~ exp(2J/T). Start- 
ing with a rectangular shape of size (L — 6) X (L+S), with 
8 oc L (5/L small), n ~ LS particles must be displaced 
from long to short edges in order to relax the deforma- 
tion. This takes a time proportional to n 2 , the motion 
being diffusive. Then we have t eq ~ n 2 T ~ exp(2 J/T)L 4 , 
namely 

= 4. (1) 

The equilibrium shape is then stabilized by a more fa- 
vorable surface to volume ratio. We have computed the 
value of (3 by means of numerical simulations of the equili- 
bration of single droplets with vacancy dynamics. Defin- 
ing R x and R y as the gyration radii along the x and 
y direction respectively, t eq is identified with the time 
needed for an initially rectangular droplet to recover the 
isotropic shape, namely R x (t eq ) / R y (t eq ) — 1. The results 
at T = 0.2 for droplets smaller than L C (T = 0.2) w 41, 
shown in Fig. ||, are consistent with Eq. ([!]). 

For rough interfaces one can instead consider an ana- 
lytic profile and invoke the classical result by Mullins that 
the equilibration time is proportional to L 4 jig]. How- 
ever, in the present context this result must be modified 
to take into account the fact that the vacancy spends 
part of the time in the bulk of domains and during such 
intervals no dynamics occurs on the interfaces. Then we 
have t eq ~ L 4 /p, where p is the fraction of time a va- 
cancy stays on domain walls. At low temperatures p « 1 
(for this reason p has not been taken into account when 
dealing with flat interfaces) whereas p ~ exp(2 J/T)I/ _1 
with rough walls. In summary we obtain 

= 5. (2) 

Numerical simulations at T = 0.6 (Fig. |]) are in good 
agreement with Eq. (||). 

Relationship between (3 and z. The vacancy- 
mediated phase-separation changes the relationship be- 
tween f3 and z with respect to Kawasaki dynamics. The 
crucial observation is that only one vacancy acts on J\f 
features of size L(t). As the typical size grows the num- 
ber of objects to be relaxed is reduced and this speeds up 
the separation process. The connection between t eq and 
the growth exponent z is different for equal and unequal 
compositions. In the symmetric case domain morphology 
is interconnected and the kinetics proceeds via the mo- 
tion of particles along the edges away from more curved 
regions. A feature of size L is relaxed to equilibrium in 
a time t eq ~ L@ , provided the vacancy acts on it. How- 
ever, since Af ~ L~ 2 such features must be equilibrated 
by a single hole, the time needed to relax them all is 
Aft e q. This relaxation causes an increase of the typical 
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size AL ~ L. Then dL(t)/dt ~ AL/(Aft eq ) ~ i 3 " and 
hence 

2 = ^-2 (3) 

Using Eqs. ([!]) and (||) we obtain z = 2 and z = 3 for low 
and intermediate temperatures, consistently with the re- 
sults reported in Fig. |l|. Notice that the exponent z = 3 
observed at T = 0.6 could be interpreted in terms of the 
usual evaporation-condensation as well. However, this 
mechanism alone is not consistent with the exponent ob- 
served at the same temperature in the quench with un- 
equal compositions. 

In the asymmetric case the minority phase forms non- 
percolating domains dispersed in a matrix of the other 
component. In this case growth proceeds via diffusion 
and coalescence of entire droplets. From the knowledge of 
(3 one can simply evaluate the diffusivity D of a domain of 
size L produced by particles moving along its edges. The 
relaxation of a droplet to equilibrium requires a number 
of molecules of order L 2 to travel diffusively a distance of 
order L; the mean square displacement x 2 of the domain 
centre of mass is then of order L. Therefore D ~ x 2 /t eq ~ 
L~ a , with 

a = (3-1. (4) 

This relation has been checked numerically by computing 
the droplet diffusivity for different L. The results, shown 
in the inset of Fig. ^, are consistent with Eq. (0) for low 
temperature. For T = 0.6 the agreement is poorer, pre- 
sumably because evaporation and condensation events 
contribute also to the motion of the droplet. The expo- 
nent a can, in turn, be related to the coarsening growth 
law by a generalization of the Binder-Stauffer argument 
fj"7f to vacancy- mediated dynamics: Two droplets sep- 
arated by a distance L take a time t coa ~ L 2 /D to 
coalesce. As above, since only one vacancy is active, 
the time needed for the coalescence of all droplets is 
Aft coa , N ~ L~ 2 being the total number of domains. 
When this happens the characteristic domain size in- 
creases of a quantity AL ~ L. Then, using Eq. (||) one 
has dL/dt ~ AL/{Mt coa ) ~ L 1 -" and 

z = a = /3-l. (5) 

With the values of Eqs. ([!]) and (0) we obtain z = 3 
and z = 4, consistently for low T with the exponents in 
Fig. [| For T = 0.6, z w 3.6 consistently with the value 
of a determined numerically at the same temperature, 
indicating that bulk diffusion weakly contributes. For 
the sake of clarity, we summarize in Table [| the exponents 
derived via scaling arguments. It must be stressed that, 
with the exception of the symmetric quench with rough 
interfaces, all other exponents are preasymtotic and will 
crossover to 1 /3 asymptotically. 



In conclusion, we have shown, by means of Monte- 
carlo simulations and scaling arguments, that the phase- 
separation of the Ising model with vacancy-mediated 
dynamics exhibits well-defined preasymptotic regimes, 
where coarsening is governed by interfacial processes. 
This results in the unexpected feature of faster growth 
at lower temperatures. In particular, for symmetric 
quenches, the effective growth exponent attains the un- 
usually large value 1/2. This whole phenomenology is no- 
tably different from what occurs in the case of Kawasaki 
dynamics. We expect these results to hold also with 
larger hole concentrations, qualitative changes occurring 
only when interfaces start to be saturated by vacancies. 
Our scaling arguments can naturally be extended to the 
three dimensional case, but compelling numerical evi- 
dence would require a huge computational effort. The 
vacancy-mediated dynamics is expected to differ from the 
Kawasaki one also in such a case. Moreover, since the 
roughening temperature Tr is finite in d = 3, the faceted 
nature of interfaces will influence even the asymptotic 
properties of the coarsening process for T < Tr. 
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Interfaces 


Symmetric quench 


Asymmetric quench 


Faceted 
Rough 


1/2 
1/3 


1/3 
1/4 



TABLE I. Value of the growth exponent 1/z for the differ- 
ent types of quench. 
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FIG. 1. Double logarithmic plot of L(t) vs t for symmetric 
quenches at various temperature. Data are averaged over 20 
realizations. 
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FIG. 2. Double logarithmic plot of L(t) vs t for quenches 
with unequal compositions at low and intermediate tempera- 
tures. Data are averaged over 20 realizations. 
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FIG. 3. Schematic representation of the fundamental pro- 
cess leading to relaxation of a faceted droplet. 
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FIG. 4. Main: Size dependence of the time t eq needed 
to equilibrate an initially deformed droplet. Best fits yield 
(3 — 3.8 and /3 = 5.1 at T — 0.2 and 0.6 respectively. Inset: 
Size dependence of the diffusivity D of a single droplet. Best 
fits yield a = 3 and a = 3.6 at T = 0.3 and 0.6 respectively. 
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